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Abstract 


The transmission of radio waves from the earth to a receiver in a 
spacecraft that is moving behind a planet makes it possible to observe 
the variations in the amplitude, phase path, or group path caused by the 
planetary ionosphere. It is shown that measurement of any one of these 
quantities may be used to determine the radial electron density distri- 
bution in that part of the planetary ionosphere which is probed by the 
signal. Measurement of two of these quantities simultaneously will 
reduce uncertainties due to changes in the interplanetary medium or 
the earth T s ionosphere during the experiment. Separate radial electron 
density distributions can be obtained from the measurements during 
immersion and emersion and these are related to the spatial properties 
of the atmosphere on the planet . It is shown that the lower atmosphere 
as well as the ionosphere of a planet can be studied by using so high 
a frequency that the signal is not influenced by the ionized component 
of the atmosphere. 

It is also shown that the statistical properties of a signal 
reflected from a planetary surface can be related to the surface statis- 
tics when the surface roughness scale is larger than the wavelength. A 
relation is found between the surface-autocorrelation function and the 
time-autocorrelation function of the reflected signal, a generalization 
of results obtained previously for monostatic radar. When a circularly 
polarized wave is transmitted, it is shown that the mean value of the 
dielectric constant and the conductivity for the planetary surface may 
be found from the shape and orientation of the polar izat 
the reflected signal. 
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I . INTRODUCTION 


Ground-based radar equipment is being used to study the moon and 
the closest planets. However, there are severe limitations associated 
with having both the transmitter and the receiver on the earth. These 
limitations can now be overcome, however, since it has become possible 
in recent years to send equipment to the moon and the closest planets. 

In the following chapters techniques are studied for exploring the 
moon and the planets using bistatic radar. Bistatic radar differs from 
monostatic radar in that the transmitter and the receiver are situated 
at different locations. 

One can transmit radio waves from the earth and receive them in a 
spacecraft moving behind a planet. If the planet is surrounded by an 
ionosphere, the signal will be influenced during the occultation. The 
changes in the signal due to the planetary ionosphere can be used to 
determine the electron density distribution and possibly some other 
interesting properties of the planetary ionosphere. The solution of 
this problem is presented in two steps. First, the changes in the 
signal are found in terms of the electron density profile of the plane- 
tary ionosphere (Chapter II) . Later (in Chapter V) the problem is re- 
versed and it is shown how the electron density distribution can be 
determined . 

More specifically, Chapter II is used to develop high-frequency 
approximations (transmitter frequency much higher than the maximum 
plasma frequency in the planetary ionosphere) for amplitude, phase, 
etc., for the waves received behind a planetary ionosphere with known 
electron density distribution. Why the high-frequency case is of such 
importance will become clear at the end of the study. 

The results derived in Chapter II are applied to the exponential 
and the Chapman ionospheric model in Chapters III and IV, respectively. 
The results are represented graphically in a normalized form so that 
they can be easily used when the exponential or the Chapman model 
applies. It is also shown that more complex ionospheres can be studied 
by considering them as being constructed from a combination of these 
basic models. 
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Chapter V is mainly devoted to the development of processes that 
can be used to determine the radial electron density distribution in 
that part of the planetary ionosphere which is probed by the signal. It 
is shown that measurement of either the amplitude, the phase path, or 
the group path during the occultation can be used to determine the 
electron density profile. 

Radio waves reflected from the planetary surface will also be 
received in the spacecraft while it is in the vicinity of the planet. 

It is, of course, necessary to be able to distinguish between the direct 
and the reflected signal. Chapter VI shows that this is possible, and 
it also shows that the reflected signal may be used to determine sta- 
tistical properties of the planetary surface, such as roughness- 
autocorrelation function, mean dielectric constant, etc. 
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II. PROPAGATION OF HIGH-FREQUENCY RADIO WAVES IN PLANETARY 
IONOSPHERES WITH KNOWN ELECTRON DENSITY DISTRIBUTION 


This chapter is devoted mainly to the study of amplitude, and of 
the phase and group paths for waves refracted in a planetary ionosphere 
with a known electron density profile. 

Later (in Chapter v) the problem is reversed in order to determine 
what can be learned about the ionosphere when the properties of the 
transmitted and the received waves are known. This is analogous to the 
circuit engineering problem of finding the two-port network equivalent 
when input and output are given. 

We will assume that a circularly polarized wave is transmitted from 
the earth and that the transmitted frequencies lie well above the maxi- 
mum plasma critical frequency, both for the earth's ionosphere and for 
the planetary ionosphere to be explored. Under these conditions the 
magnetoionic theory shows that the effect of magnetic fields on ampli- 
tude and phase can be neglected. 

The results derived for phase path and amplitude are also valid 
when a linearly polarized wave is transmitted from the earth. In the 
latter case one can also observe Faraday rotation. Assuming that the 
Faraday rotation due to variations in the earth 1 s ionosphere can be 
neglected or subtracted out, transmitting linearly polarized waves 
yields a means of exploring the planetary magnetic field (if there is 
any). This problem is discussed later. 

The geometry is illustrated in Fig. 1. That part of the trajectory 
along which occultation can be observed may often be well approximated 
with a straight line. The y axis goes through the center of the planet 
and the transmitter on the earth. The x-y plane is chosen parallel to 
the fly-by trajectory. 

For a spherically symmetric ionosphere, it is convenient to rotate 
the points along the trajectory around the y axis until they end up in 
the x-y plane. In this way a new curve, which here will be called the 
equivalent trajectory, is obtained in the x-y plane. In the spherically 
symmetric case, it is sufficient to study the phase and amplitude of the 
waves in the x-y plane, which will then also give phase path and ampli- 
tude along the true trajectory. 
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FIG. I. GEOMETRY DURING THE OCCULT AT ION. 


In the following sections the ray theory is applied to analyze the 
propagation through the planetary ionosphere. The refracting medium is 
inhomogeneous but isotropic. The energy is therefore propagating along 
the wave normal and this direction can be found from Snell* s Law. 

Necessary corrections to the ray-theory solution are taken up later. 

A. RAYPATHS 

Refraction in spherically symmetric ionospheres is discussed first 
because the effect of nonspherical perturbations in the electron density 
distribution is small. A polar coordinate system can be used to describe 
the raypaths in the x-y plane when the ionosphere is assumed spherically 
symmetric . 
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The differential equation for the raypaths can be derived in the 
following way. From the triangle ABC in Fig. 2 it is seen that 


r 


dr 


tan i 


( 2 . 1 ) 


To obtain an explicit expression for the rays in terms of r and 
a relation between i and the refractive index p(r) is needed. It is 
seen from Fig. 3 that for an ionosphere built up of shells (spherical 
layers) with constant p, one has 


r i Sln i l = r 2 Sin ?2 
Sin ?2 = ^2 Sln i 2 
r 2 Sin i 2 = r 3 Sin ?3 


(law of sines) 
( Snell's Law) 
(law of sines) 


p „ sin y = Li sin i (Snell's Law) 

n-1 'n n n v ' 

The product of these equations gives 

p r sin i = p r sin i 
11 1 n n n 

or when (a(r) is a continuous function of r, 

pr sin i = p r = a ( 2 - 2 ) 


where r^ is the radius of closest approach for the raypath and p^ is 
the refractive index at this point. The distance between the center of 
the planet and the raypath asymptote is denoted by a. Equation (2.2) 
is sometimes called Bouguer's rule. 

Combining Eqs. (2.l) and (2.2) yields the differential equation for 
the raypath 


dr 




/ ri r\ \ 


This is a fairly well known result [Ref. 1]. 
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Assuming that the effects of the magnetic field and the collisions 
are negligible, we get from the magnetoionic theory 


|i 2 ( r ) = 1 - X(r) 


(2.4) 


where 


X(r)=^N(r) (2.5) 

r 

Here, f is the frequency and N(r) is the radial distribution of 
electron density in the planetary ionosphere. 

By combining Eq . (2.3) and Eq. (2.4) and integrating, one obtains 


4 ) 



i - x(iO 
1 - 



( 2 . 6 ) 


where <J> = o for r = r . 

o 

Assuming that the electron density profile is known, one can use 
Eq. (2.6) to determine raypaths for different r^. This family of ray- 
paths (Fig. 4) may be used to find the ray pictures when some radiating 
source is present inside or outside the ionosphere. For instance, when 
the source is far away from the planet, all the rays should be rotated 
about the origin so that one of their asymptotes points towards the 
source . 

In most fly-by experiments, it is expected that the trajectory will 
lie well outside the ionosphere, and this makes it possible to simplify 
the procedure. Instead of finding the total raypath, it is sufficient 
to determine the asymptotes. The geometry is shown in Fig. 5 where: 



(2.7) 


and 


^ cos (2\|r + g ) R p 

The new geometrical quantities involved are shown in Fig. 5. 


( 2 . 8 ) 
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FIG. 


N 


PLANETARY SURFACE 


FIG. 4. FAMILY OF RAYPATHS . 


y 



5. IMPORTANT RAYPATH PARAMETERS. 
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We will now introduce a new variable OL given by 


r / 


(2.9) 


a - arccos 

Now one obtains for the total angle (2\[f) that the rays are refracted 

it/2 


2<|r = 2 


f 


r 

2 

/I - 

1 

1 

o 

j 

1/2 

\ da 

(2.10) 

< 

2 

x(r ) cos a - X( r /cos a) 
li° ° 



.2 

_ sin a 

J 



This equation is correct for a wide range of frequencies as long as 
collisions and magnetic fields are absent. However, we are mainly 
interested in the high-frequency case (x« l), and it would be 
desirable to simplify Eq . (2.10). Looking at 


5(a,r ) = X(r ) 
o' o 


2 

cos a 


X(r /cos a) 
o 7 

«‘-o> 


. 2 

sm a 


( 2 . 11 ) 


one sees that | 5| is less than 1 for a greater than 0 because 

X(r ) « 1. For I a| « 1 we have 

o .11 


X 



X(r o ) ♦ f a 2 x'(r ) * 


Thus for small a: 


S(a, r o ) = - x(r Q ) 


1 


r o X ’ (r o ) ' 
* 2 


which also turns out to be small compared to 1 for the ionospheric models 
we will consider. Using this, we can expand the integrand in Eq . (2 .10) 
and, maintaining only the first-order terms in X, we obtain: 

x(r /cos a) 
o 

X(r ) 

2 * 1 1 = x(r Q ) / 2“^ da (2.12) 

O sin a 


n/2 1 


when |&(a,r o )| « 1. 
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For the cases we shall consider, 2\|r turns out to be less than 1 or 2 
deg, and Eq . (2.12) usually provides good accuracy. 

Equation (2.12) shows that there is a linear relation between 2 \|r 
and X. If we have an ionospheric model [x(r Q /cos a)] built up of 
several layers [ X^ (r^/cos Of)], this can be written 

x p_) X (—2—) (2.13) 

\cos OCJ / J v ^ cos Oi J 

Applying Eq . (2.12), we find for the refracted angle: 

2 \| r 2 t v (2.14) 

v 

where 2\|f is caused by the v** 1 layer, etc. In other words, super- 
position applies, and refraction due to each layer may be treated as if 
it were the only layer present . 

In Chapters III and IV, exponential and Chapman-distributed electron 
density profiles are assumed. More complex ionospheres can be studied 
by considering them as being constructed from a combination of these 
basic models. 

We have so far assumed a spherically symmetric ionosphere and, on 
this basis, determined the raypath. However, one can also start by 
assuming an ionosphere with nonspherical perturbations. Let us look at 
the wavefronts as a wave propagates through a planetary ionosphere for 
the high-frequency case illustrated in Fig. 6. The wavefronts become 
slightly perturbed as the wave propagates through the ionosphere. The 
phase at point P in Fig. 6 is proportional to the integral of p along 
the raypath. For the case shown, it does not matter whether one inte- 
grates p along the raypath or along the straight line H . (This is 
only true when the point P is close behind the ionosphere and when the 
refraction is small.) Using this argument, we see that /( 1 - p) dy 
taken along l gives the perturbation on the wavefront in the y direction. 
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FIG. 6. PERTURBATION OF THE WAVEFRONTS. 


The angle the wave normal makes with the y axis ( 2\|r ) is equal to 


minus the slope of the wavefront, and we therefore obtain: 


through 


medium 


(2.15 


where the gradient in [i is assumed negligible perpendicular to the 
x-y plane. 

Equation (2.15) states that the angular deviation of the rays in the 
high-frequency case is equal to the gradient in the integrated phase 
refractive index, \i . This result has, for instance, been used to 
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determine horizontal gradients in the total ionization of the earth r s 
ionosphere by measuring the ionospheric refraction of waves from radio 
stars [ Ref . 2 ] . 

Let us apply Eq . (2.15) to the spherically symmetric case and 
compare this with the results we obtained previously. Equation (2.4) 
gives in the high-frequency case: 

H(r) = 1 - | X(r) (2.16) 


thus 



dy 



dy 


= - 1 h / x(r) dy 

- CO 


After introducing x 


r and y = r tan a, we obtain: 
o J o ' 



o 




x(r /cos a) 
o' 

cos^a 


da 




r da 
o 

cos^a 


The last integral can be integrated by parts giving 



X(r /cos a) 
o' 

cos 2 a 



1 

sin^a 


da 


X( 


r ) lim ( — : — 

° a-o \ sin a 


a cos a 


Further manipulations give 

x(r^) " X(r^/cos a) 


2 * -f 


^s in^a 


da + X(r ) lim 
o' 

a 


o( 


sin a cos a 


cot a 
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where the limit is zero. By comparing this expression with Eq . (2.12), 
we see that the two approaches yield the same result, and we have there- 
fore justified the simple arguments leading to Eq . (2.15). 

We have so far limited the discussion to refraction caused by 
gradients in the ionization. However, it is worthwhile noting that the 
same results can be used to treat refraction in the lower atmosphere. 

For the refraction index in this case, we can use [Ref. 3]: 


M' 


1 + 10 



380, 000e\ 
T 2 ) 


(2.17) 


where P is the partial pressure of dry air in millibars, e is the 
partial pressure of water vapor in millibars, and T is the tempera- 
ture in degrees Kelvin. This corresponds to replacing Eq . (2.5) by 


X = 


2-10 


-6 



380 , 000e\ 
T 2 ) 


(2.18) 


When T is constant throughout the atmosphere, the pressure can be 
expressed as an exponential function of height. We can therefore use 
Chapter III to take into account the influence of the atmosphere. 


B. PHASE PATH (s ) AND GROUP PATH (s ) 

P g 

The ray picture is useful in the visualization of the propagation, 
but we are more interested in those quantities which can be measured, 
such as phase and group paths. The phase path of the signal is found 
by integrating the phase refractive index p along the raypath from the 
transmitter (or some reference surface), through the planetary ionosphere 
and to the receiver. 

There are two reasons why the phase path changes as the transmitter 
is occulted by the ionosphere: 

1 . The electrons advance the phase and therefore tend to make the 
phase path shorter. 

2. The raypaths are curved due to refraction and this tends to 
increase the phase path. 

It is necessary to study both these competing effects in order to draw 
any conclusions about how the phase path is changing. 
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Previously we indicated that the phase path could be found by inte- 
grating \i along the raypath. This method is sufficient when the 
propagation takes place over relatively short distances, such as between 
two points on the earth T s surface. However, it is not a practical 
approach in our case because for most fly-by trajectories, one has to 
integrate over such huge distances that the effect of the planetary 
ionosphere more or less disappears. One can improve the accuracy, when 
a digital computer is available, by using double-precision arithmetic. 
However, since the use of double-precision arithmetic slows down the 
execution of the program considerably, it is not a very attractive 
solution from an economical point of view. Hence, it is desirable to 
look for a more sophisticated method for the calculation of the phase 
path . 

Our primary interest is not so much the total phase path itself, 
but rather the changes in the phase path due to the planetary ionosphere. 
These changes would be fairly simple to compute if the raypaths were the 
same with and without the ionosphere. In the latter case, it would only 
be necessary to integrate the electron density along straight lines. 
However, we must look deeper into this because the rays are curved. 

Defining s as the decrease in the phase path due to the presence 
P 

of the planetary ionosphere, we can write: 



along along 

straight raypath 

line 


where the second integral is taken along the actual raypath, and the 
first integral is taken along the raypath connecting transmitter and 
receiver when the planetary ionosphere is (imagined) removed. The 
parameter c is the free-space phase velocity and f is the frequency. 
Thus, 3^ gives the phase difference in cycles. 
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In the high-frequency case, Eq. (2.19) gives 


§( / ds - / ds + 5 / Xds )= S P 

along along along 

straight raypath raypath 

line 


+ s r 


(2.20) 


where 



raypath straight 

line 

s 2 = h J * ds 

along 

raypath 


(2.21) 


(2.22) 


The group path is of interest when one transmits pulses instead of 
cw. For the increase in group path (s^) °* ue t ^ le presence of the 
planetary ionosphere, one obtains 


s 

g 


S 1 + S 2 


(2.23) 


when magnetic fields and collisions can be neglected. The group time 
delay due to the planetary ionosphere is (s^/f) with this notation. 

Equation (2.2l) gives s^ in terms of two integrals; the magnitude 
of these integrals is very large compared to their difference. It is 
therefore not practical, as described earlier, to compute the two inte- 
grals and then subtract them. Instead, one can compute s^ in two 
steps. The first step consists of finding s as if the propagation 
was along the raypath asymptotes. The actual raypath between the trans- 
mitter and the receiver may be shorter or longer than the distance 
measured along the raypath asymptotes, and it is the purpose of the 
second step to correct for this difference (As n ) . 

X 

One obtains, when the rays can be considered parallel before they 
enter the planetary ionosphere: 

s l = I ( x . a ) tan * + As^ (2.24) 
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where x designates the x-coordinate of the receiver for the equivalent 
trajectory, and a and t are given by Eq . (2.2) and Eq . (2.12), 
respect ively . 

Computation of As^ requires integration. Use of a polar coordinate 
system to find As^ is not practical because the errors get larger as 
the distance to the planetary center is increased and the process does, 
for this reason, not converge. Figure 7 shows a more convenient coordi- 
nate system giving fast convergence when £ is used as a dummy variable 
of integration. It turns out that As^ represents a higher order 
correction to s^ and can often be neglected, as later results will show. 

The phase advance or group delay due to the electrons is given by 
Eq. (2.22). Again using the polar coordinate system and a as a dummy 
variable, we obtain: 



FIG. 7. A CONVENIENT COORDINATE SYSTEM 
FOR CALCULATION OF As x . 
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= 


X(r /cos a) [ 

•J 


cos 




sin a + 


2 

cos a 


(**) 

\ da/ 


,i/ 2 


da 


(2.25) 


where 


= [1 - X(r )] 
(da/ v o /J 


1 + 


X(r^) cos a - X(r Q /cos a) 


i-l 


si 


in 2 a 


(2.26) 


Expanding the integrand in Eq . (2.25) and maintaining only the first- 
order term in X gives 


s 


2 



x( r^/cos a) 
cos^a 


da 


(2.27) 


which is the same as integrating X along a straight line through the 

ionosphere. It is important to notice here that while refraction can be 

neglected as far as computation of s goes, one cannot at all neglect 

the influence of refraction on s or s when the receiver is well 

P g 

behind the ionosphere. 

Combining Eq . (2.15) and Eq . (2.27), one obtains: 


2\|f = - 


o dS 2 (r o ) 


f dr 

o 


(2.28) 


which will prove useful in later chapters, 


The time rate of change in s^ gives the doppler shift due to the 

ionosphere. One may measure ds /dt in the following way. Assume that 

P 

two harmonically related frequencies f^ and f ^ are transmitted from 

the earth. Thus f m = f , where m is an integer larger than 1. 

i 

Before occult at ion takes place, one will receive the frequencies 

f^[l - (v^/c)] and f ^ [l - (v^/c)], where v^ is the velocity at which 

the transmitter and receiver move apart . The lowest frequency will in 

addition be shifted by Af . during occultation. We will assume that 

1 


f 2 is so high that this wave propagates undisturbed through the 
planetary ionosphere. After mixing of f A 1 - (v /c)] and the m L 


th 


mAf . 


This 


harmonic of f,[l - (v /c)l - Af., one can filter out 
1 v r 1 1 

yields s^ after integration. It is here assumed that the effect of 
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changes in the earth’s ionosphere and in the interplanetary medium can 
be neglected or predicted. 

The total electron content of the earth’s ionosphere is about 

17 2 

2*10 elect rons^meter in the daytime. Ionospheric irregularities 
cause fluctuations of about ±2 percent in the daytime content [Ref. 4], 
This variation corresponds to a change in of about ±10 wavelengths 

at frequency 50 Me. Nighttime variations in s^ would be smaller. 

Even less is known about the density fluctuations one may encounter 
in the interplanetary medium. However, the Stanford Center for Radar 
Astronomy has radio propagation experiments to be included on Pioneer I 
and II which will help answer these questions. 

Group-path measurements may also be used. This can be done by trans- 
mitting two pulse trains of frequency f and f 0 , respectively, and 
observing the change in the time delay between pulses of different 
frequency . 

Other possible measurement techniques are discussed by Eshleman, 
Gallagher, and Barthle [Ref. 5], 

C. REFRACTION GAIN (G ) 

r ' 

Another measurable quantity of great interest is the amplitude of 
the signal. Rather than compute the absolute value of the amplitude, 
we will find the changes in the amplitude due to the presence of the 
planetary ionosphere. 

Let us define the "refraction gain" ( G r ) as 

S 2 

G = 10 log — (db) (2.29) 

r S 1 

where is the Poynting vector one would have at the position of the 

receiver without any planetary ionosphere, and S is the Poynting 
vector at the receiver in the presence of the ionosphere. Figure 8 
illustrates, considerably exaggerated, the effect of refraction on the 
amplitude. We will assume that the gain of the transmitting antenna can 
be considered constant over a solid angle, including that part of the 
fly-by trajectory along which the occultation is observed. This assump- 
tion implies that, in the absence of any planetary ionosphere, the power 
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2 



FIG. 8. REFRACTION GAIN DUE TO PLANETARY 
ATMOSPHERE. 
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density is the same at positions 1 and 2 in Fig. 8. Hence one can relate 

S and S to geometric quantities on the ray picture 
1 2 

Vl^l = S 2 r 2 Ar 2 


which gives 


G 

r 


10 log 


r i Ar i \ 

r 2^2 ) 


(2.30) 


We can simplify Eq . (2.30) further because the angle by which the rays 
are refracted (2\|/) is very small in the high-frequency case. Let the 
coordinates of the receiver be (x,y). We have, using quantities defined 
in Fig. 5: 

r 2 = R p + q + (^ + P) y 

and 


r l cos (3 


+ Py 


After introducing r^ = x, one obtains 



(2.31) 


where only first-order terms in 24 f and p are maintained. We can 
also write 


_ A d2t A Dy 

Ar 0 = Ar + - — Ar_ — 

2 1 dr 1 D + y 

o 


Combining Eqs . (2.30), (2.3l), and (2.32), one obtains 


G - 10 log 
r 


1 - 2t - 

x 


1 + 


d2t Dy 


dr D + y 
o 


Close behind the ionosphere one has 


2\|r ~ « 1 

x 


d2^ Dy 
dr^ D + y 


« 1 


(2.32) 


(2.33) 
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and this gives for the gain 


G « 4.34 
r 



(2.34) 


where one can often neglect the term (2i|ry/x). 

We have here based the definition of the gain on the ray theory. 

This implies that does not describe the amplitude variations close 

to shadow boundaries or caustics. Approximations for the amplitude can 
in these regions be found from the wave theory as shown in Appendix A. 


D. CAUSTIC FORMATION 

The distance (y c ) behind the ionosphere at which the rays cross 
each other can be found by setting the denominator of Eq . (2.33) equal 
to zero: 


■m Dy c 

dr D + y 
o 


0 


or 

„ m (.m. iV 1 
c l dr o D ) 


Caustics are formed only when 


(2.35) 


d2\|f 1 

dr D 

o 


(2.36) 


for some finite interval in r , while d2\|f/dr > - 1 /D tends to 

o ' o 

defocus the energy. 
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III. APPLICATION TO AN IONOSPHERE WITH EXPONENTIAL ELECTRON DENSITY 
DISTRIBUTION 


One can, for planetary atmospheres in hydrostatic equilibrium, show 
that the electron density will be an exponential function of height when 
the density in the atmosphere is sufficiently low. This is a special 
case of the Chapman theory [Ref. 6]. 

We will in this chapter use 


x( r ) = X^ exp 



(3.1) 


where X is the normalized electron density at the surface, H is the 
s J 

scale height in the atmosphere, and R^ is the planetary radius. 

The complete Chapman formula will be used in the next chapter. 


A . RAYPATHS 

The raypaths can be determined by combining Eq. (2.6) and Eq . (3.l). 
However, it is sufficient for our purpose to know the raypath asymptotes 
and these are determined by Eqs . (2.8) and (2.12). Using Eqs . (2.12) 
and (3.l), we obtain the angle of refraction (2\(r): 


2 'if i 

r / 1 - (' - cos 0)] 

da 

(3.2) 

J 

sin^a 

O 




This is shown in Fig. 

9 as a function of r /2H. 

o f 




Equation (3.2) gives the high-frequency approximation for 2^ . 
However, we would also like to know what accuracy we can expect from the 
high-frequency approximation. This can be seen by comparing it with the 
exact solution. Combining Eq . (2.10) and Eq . (3.l), we find the exact 
solution: 


24 / 2 



o 


1 - x < r o) 


1/2 


cos 2 « - exp [- h; (1 - ^h^)]| 


da 


1 . x(r o ) 


sin^a 


( 3 . 3 ) 
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10.0 



FIG. 9. NORMALIZED ANGLE OF REFRACTION FOR THE 
EXPONENTIAL IONOSPHERE. 


The right-hand side is in this case a function both of r Q / 2H and 
X(r Q ). It is therefore necessary to choose some particular value of 
X(r ) in order to compare the two solutions given by Eq . (3.2) and 
Eq , (3.3) respectively. With x(r^) ~ 0*01 in Eq . (3.3) one obtains 
the lower curve in Fig. 9. 

Figure 9 may give the false impression that the accuracy in 2^r is 
decreasing with increasing r^ when the high-frequency approximation is 
used. However, this is not true because X(r^) falls off very fast with 
increasing r . 
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As an example, let us use H = 50 km and X =0.01 for the moon. 

s 

This gives r Q /2H ^ 17 for the grazing ray. From Fig. 9 we find that 
this corresponds to an error in 2 \|r of about 3 percent . The error is 
less for rays missing the surface. 

It may seem surprising that the curves on Fig. 9 go through the 
origin. This may occur in two different situations: 

1. In the r =0 case, there is no planet in the center of the 
ionosphere. The ray passing through the center of the ionized 
blob is not bent in any direction. 

2. In the l/H = 0 case, 2 \|r is also zero because there is no 
gradient in the ionization. 

The curves in Fig. 9 have been found by integrating Eq . (3.2) and 
Eq. (3.3) on a digital computer. Analytic techniques can also be 
employed. One can, for instance, expand the integrand into the follow- 
ing form: 


exp (- ka 2 ) ^ ^ c^ a v (3.4) 

v^or^e, ... v 

This form enables us to evaluate the integral when the upper limit ( Jt/2 ) 
can be replaced by infinity. For 2^ one can derive 



(3.5) 


A numeric example has been worked out in order to illustrate the 

theory developed so far. The example assumes an exponential model for 

the lunar ionosphere with R = 1738 km, H = 50 km, f = 50 Me, and an 

P 3 

an electron density at the surface equal to 200 electrons/cm 

(X = 6.44*10 6 ) . The corresponding ray picture is shown in Fig. 10. 

The scale is different along the two axes, and the moon will therefore 

look like an ellipsoid in this space. The waves are assumed to be 

transmitted from the earth. 


B. PHASE PATH (s ) AND GROUP PATH (s ) 

P g 

It was shown in Sec. IIB that the changes in the group path and the 
phase path, due to the planetary ionosphere, can conveniently be ex- 
pressed as the sum of, and the difference between, s 2 and s^ 
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KILOMETERS 



TRANSMITTER (ON EARTH) 


L 

o 


J L 

100 200 


_L 

300 


400 


KILOMETERS 


FIG. 10. REFRACTION IN A LUNAR IONOSPHERE AT 50 Me WITH THE TRANSMITTER 

ON THE EARTH. (Assumed electron density: N(h) = N exp (-h/2H), 

ft / 2 < max x ' 7 

whprp N = 2*10° electrons/m ; H = 50 km./ 
max 
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s 2 is approximately proportional to the electron 


s^ is proportional to the 


respectively. Here 

density integrated along straight lines, and 
difference in path length with and without the planetary ionosphere. 

The high-frequency approximation for s 2 is given by Eq. (2.27). 
Combining this with Eq. (3.l), one obtains: 

.*/ 2 exp f ~ / 1 - — i — )] 

1 2H \ cos a / J 


CS 2 f 

fr X(r ) / 

o ' o' J 


cos 


w 


da 


(3.6) 


The result is shown in Fig. 11 as a function of r Q /2H (upper curve). 

Equation (3.6) gives the high-frequency approximation for s^. We 

also like to know what accuracy we can expect from the high-frequency 

approximation. One can determine the accuracy of Eq. (3.6) by comparing 

it with the exact solution given by Eq. (2.25) and Eq . (2.26). The 

exact solution can be normalized in the same way as Eq. (3.6). However, 

the right-hand side will be a function both of r /2H and x(r ). It 

o N o' 



FIG. 11. NORMALIZED s FOR THE EXPONENTIAL IONOSPHERE. 
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is therefore necessary to keep x(r^) constant in order to compare the 
exact solution with the high-frequency approximation. Choosing 
x(r^) = 0.01, one finds the lower curve in Fig. 11. 

The pole at the origin in Fig. 11 is caused by the normalization 
and does not mean that s is infinite for a ray passing through the 
center of a blob of exponentially distributed ionization. 

The curves on Fig. 11 have been obtained by evaluating the integrals 
on a computer. Analytic techniques can also be employed, as already 
mentioned in Sec. IIIA. For instance, one can show that: 


- s/flHr” X(r ) 
c o v o' 


(3.7) 


For certain intervals in i~ o /2H, one can find better approximations 
for s 2 . 

We also need to find s^ given by Eq. (2.2l). The rays can first 
be replaced by their asymptotes as described in Sec. IIB. This computa- 
tion does not require any new integration when is known. Equation 

(2.24) gives s^ when the waves can be considered plane before they 
are refracted in the planetary ionosphere. The term As^ is the 
difference between the phase path measured along the raypath and the 
phase path measured along the raypath asymptotes. 

The term As^ can be approximated in different ways. For instance, 
one can replace that part of the raypath that is closest to the plane- 
tary ionosphere by an arc of a circle. This raypath shape gives 


As., = 2 — (\|r - tan \|r ) r (l - — sec \|r ) 
1 c ' o\ r / 


sec - 1 


- | - r tx(r ) 

3c o v o' 


(3.8) 


which is usually negligible at high frequencies. 

The raypath is not really an arc of a circle close to the surface, 
and it is therefore necessary to check how much the actual shape of the 
raypath will influence the form of the expression for As^. For in- 
stance, one can replace the arc of the circle by a piece of a parabola 
with maximum curvature at the point of closest approach. (The reason 
for choosing a parabola is, of course, that the line integral of this 
conic section can be found in terms of elementary functions.) By com- 
paring the two results one finds that the difference in As^ is negligible. 
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The reason why can be approximated in such a simple way is, of 

course, that the raypaths are almost straight lines. 


C. REFRACTION GAIN (G ) 

V r 

The refraction gain (g^) was defined in Sec. IIC as the change 
in signal strength due to the refraction in the planetary ionosphere. 
Equation (2.33) shows that this gain depends on which part of the 
ionosphere the signal is propagating through and how far the receiver 
is behind the planetary ionosphere. The numerator of Eq. (2.33) 
contains 2^ which can be found from Eq. (3.5), and the denominator 
of Eq. (2.33) can be approximated in the following way: 


1 1 d(2 ^ Dy ~ 1 Dy s /p 

D + y 2 l H 


dr D + y 
o J 


l/2 

X_ /*R_\ /R - r 

/ P o 
exp' — — — 


v 


2H 


, r - R 
1 1 o p 

2H 2R + 4HR 

P P J 


(3.9) 


where 2i|r has been approximated by 

l/2 


x /hr 

2 * “ “f (“h 2 


1 r - R 

1 + 1 -JL- E 

2 R 


exp 


R - r 
_P o 

2H 


D. CAUSTIC FORMATION 

An approximation for the caustic can be found by combining Eq. (2.35) 
and Eq. (3.5). This corresponds to setting the right-hand side of Eq. 
(3.9) equal to zero. One then obtains for the y-coordinate (y Q ) of 
the caustic: 



1 

(3.10) 


This expression can be simplified further when 

« 2R . If one also approximates x with ( 
P c 


H « R and 
P 

r + y 2\|r ) , one 
o J c f 


r - R 


o p 

obtains 


where 





designates the coordinates of the caustic. 


(3.11) 
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A ray picture was shown in Fig. 10, but it did not have any caustic, 
the reason being that the curvature of the wavefronts was too large 
before the waves were refracted in the lunar ionosphere. However, the 
result will look different if the wavefronts are assumed plane before 
they are refracted in the moon's ionosphere. This is illustrated in 
Fig. 12. The ionospheric model and the frequency are the same as in the 
previous example shown in Fig. 10, but the source is at infinity. 

The caustic in Fig. 12 is given by Eq . (3.1l) with 1 /D = 0, which 
gives 


y 


c 





(3.12) 


In order to have a caustic, it is necessary that y c in Eq . (3.10) 

be positive for some finite interval in r . This requires that 

o 


X 

s 




1 

H 



(3.13) 


We note that at sufficiently low frequencies there will exist a caustic, 
but as the frequency is increased, ceases to satisfy the above 

relation, and the caustic disappears. 
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KILOMETERS 



FIG. 12. REFRACTION IN A LUNAR IONOSPHERE AT 50 Me WITH THE SOURCE AT 
INFINITY. (Assumed electron density: N(h) = N exp (-h/2H), where 

Q ry III 8.X 

N = 2*10 eleetrons/m ; H = 50 km.) 
max ' ’ ' 
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IV. APPLICATION TO AN IONOSPHERE WITH CHAPMAN-DISTRIBUTED ELECTRON 
DENSITY 


The Chapman theory [Ref. 6] gives the electron density as a function 
of solar zenith angle (x) and radial distance from the center of the 
planet (r). Normalizing the Chapman electron density distribution 
using Eq . (2.5), one obtains: 


\cos a max 


exp 


1 + 


R +h 
_P o 

H 


H cos CL 


- sec 


X expl 


/R +h 
1 p o 


H 


H cos QL 


(4.1) 

where r^/cos a represents the distance to the center of the planet (r) 
[using Eq.(2.9)], X is the solar zenith angle, H is the scale height 


in the atmosphere, R is the radius of the planet, 


is the height 


of the maximum electron density at X = 0, and X is the maximum 

in the normalized electron density. The Chapman ionosphere is not 
spherically symmetric. 

The signal transmitted from the earth and received in a space probe 
behind the Chapman ionosphere will only propagate through part of the 
planetary ionosphere. The variations in the electron density with 
longitude and latitude can be neglected as a first approximation in that 
part of the ionosphere which is probed by the waves. This corresponds 
to setting X in Eq. (4.l) equal to a constant X Q * For X Q we will 
use the solar zenith angle at the point where the ray passing through 
the receiver is closest to the planet . 

Equation (4.l) can now be rewritten in a more convenient form: 


( — ) = X \ / cos X * exp 

\cos CL/ max ^o 


1 + K - 


H cos a 


- exp ( K 


H cos a 

(4.2) 


where 


R + h 
p o 

H 


+ In sec X 


t . _ \ 


Equation (4.2) represents the normalized electron density distribu- 
tion in a spherically symmetric ionosphere. The maximum X is 
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X \ I cos X at a height h . The height h is given by 

max o max max 

h = h + H In sec Y (4.4) 

max o^o v 7 

It is important to notice here that we are not replacing the entire 
Chapman ionosphere with a spherically symmetric ionosphere. We are 
merely neglecting the nonspherical variations in the relatively small 
part of the ionosphere which is probed by the signal. This simplifica- 
tion will be further justified later. 

A. RAYPATHS 

Section IIA describes how one can determine the raypath asymptotes 
when the ionosphere is given. The determination of the raypath 
asymptotes mainly amounts to calculation of the angle (2\|r) that the 
rays are refracted. 

Combining Eqs . (2.12) and (4.2), one obtains the high-frequency 
approximation for 2\|r : 



The right-hand side can be considered a function of K and (h-h )/2H. 

Here h denotes the height at which the ray is missing the surface, 

thus h is equal to (r - R ). The result of integrating Eq . (4.5) is 

op 

shown in Fig. 13. 

Figure 14 gives an indication of what accuracy one can expect from 

the high-frequency approximation by comparing it with the exact solution. 

The exact solution is found by combining Eq . (2.10) and Eq . (4.2). 

Normalizing the exact solution in the same way as in Eq . (4.5), 

one obtains a function of K, (h-h )/2H, and X \l cos Y on the 

max' max ^o 

right-hand side of the equal sign. The curve representing the exact 

solution in Fig. 14 is drawn for K = 30 and X \ cos Y = 0.01. 

max o 


exp 


II 11 


+ K - 


H c os a 
in 2 a 


- exp IK 


H cos OU 


si 


da 


(4.5) 
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FIG. 13. NORMALIZED ANGLE OF REFRACTION FOR THE CHAPMAN IONOSPHERE 




[Sg gj] 


IONOSPHERE . 






B. PHASE PATH (s ) AND GROUP PATH (s ) 

P g 

It was shown in Sec. I IB how the changes in group path or phase path, 
due to the planetary ionosphere, can be expressed respectively as the 
sum of or the difference between s and s . 

Sections IIB and IIIB showed that s^ can be approximated as the 
phase path along the raypath asymptotes minus the phase path along a 


straight line connecting transmitter and receiver. Thus s^ can be 
found without integration when the raypath asymptotes are determined. 
Equation (2.27) gives the high-frequency approximation for s . 

Ci 

Combining Eq . (2.27) and Eq . (4.2), one obtains: 


cs r 


n/2 


fX r 
max o 


4 


sec X 


>/ ^11 


exp o 1 + K - 


H cos a 


exp IK 


(k - r ° ) 

\ H cos a A 


da 

cos^a 


(4.6) 

The result is shown in Fig. 15 where h is the raypath miss distance and 

h is the height of the maximum electron density, 

max 

Figure 16 shows a comparison between the high-frequency approximation 

given by Eq . (4.6) and the exact solution given by Eqs . (2.25), (2.26), 

and (4.2) for K = 30 and X Jcos X = 0.01. 

max o 

For large values of (h-h ) /2H. one can use the results in 

max ' 

Chapter III. 

The integrals in this chapter have mainly been evaluated on a digital 
computer, but analytic techniques can also be employed. The integrand 
in Eq . (4.6) can for instance be expanded as described in Chapter III. 
However, different intervals in a require different expansions, and 
it is therefore difficult to obtain good accuracy when adding and sub- 
tracting the different areas. For rays that do not pass inside the 
maximum in the electron density, one can show that 


f 

- X 
c max 


(rtHr cos 
o 


vl/2 [1/. 

>> eXP t2l" 


r o\l 

¥7] 


3 

8 


2 

3 


1 

2 


!°\1 

H /J 


(4.7) 
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Similar approximations have been found for other intervals in (h-h ) /2H, 

max 7 7 

but the accuracy of these approximations is not very good and they are 
therefore not given here. 



FIG. 15. NORMALIZED s 2 FOR THE CHAPMAN IONOSPHERE. 
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16. COMPARISON BETWEEN EXACT AND APPROXIMATE 
FOR THE CHAPMAN IONOSPHERE. 
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C. CAUSTICS 


* 


The refraction gain (g^) behind a Chapman ionosphere can be 
computed using Sec. IIC and Sec. IVA. The gain is infinite along 

caustics and this fact is made use of in Sec. IID to obtain the equation 
for the caustic [ Eq . (2.35)]. 

Combining Eq . (4.5) and Eq . (2.35), one obtains for the caustic 
ordinate behind a Chapman ionosphere: 



(4.8) 

The result is shown in Fig. 17 as a function of (h-h )/2H, where 

h is the raypath miss distance and h is the height at which the 

J max 

electron density has its maximum. We see from Fig. 17 that two 
caustics may exist behind a Chapman ionosphere. 



FIG. 17. NORMALIZED DISTANCE TO THE 
CAUSTIC FOR THE CHAPMAN IONOSPHERE. 
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D. NUMERIC EXAMPLE USING A CHAPMAN MODEL FOR THE MARTIAN IONOSPHERE 


The previous sections have shown how the refractive properties of 
the Chapman layer can be represented graphically in a normalized form. 

In this section a numeric example is given assuming that Mars has a 
Chapman ionosphere. The signals are transmitted from the earth and 
received in a spacecraft that is being occulted by the Martian iono- 
sphere. We will compute phase path and refraction gain of the signal 
received in the spacecraft. The geometry is illustrated in Fig. 18. 

Two receiver trajectories are shown. 

The radial distribution of the electron density along the x axis is 
shown in Fig. 19. We will first neglect the nonspherical variations in 
the electron density distribution, thus 

X - x c 

where Y is the angle between the x axis and the direction to the sun. 
o 

This approximation will be justified later. 

Figure 20 shows the ray picture in the x-y plane for f = 50 Me . 

The curvature of the wavefront is much smaller before it reaches Mars 
than after it has been refracted in the Martian ionosphere. The waves 
have therefore, for simplicity, been assumed plane before they are 
refracted in the Martian ionosphere. 

The atmosphere on Mars also causes refraction. Chapter III, with 
X defined by Eq . (2.18), is well suited for making estimates of the 
atmospheric effects. Reference 7 gives the necessary atmospheric 
parameters for Mars. Using these parameters one finds that the effect 
of the atmosphere can be neglected. 

Figures 21 and 22 show how amplitude and phase change when occult a- 
tion takes place along trajectory 1. The doppler shift of the signal, 
caused by the Martian ionosphere, is given by 




ds 

n 

£1 

dt 




cps 


The first factor (ds^/dx) 

distance (x - R ) in Fig. 

p 


is shown as a function of receiver miss 

3 

23. For -dx/dt = 10 m/sec we see that 
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FIG. 20. REFRACTION IN A MARTIAN IONOSPHERE AT 50 Me. (Assumed electron density: Chapman layer, 

H = 130 km, h = 320 km, N = 10 11 electrons/m 3 , sec Y = 1.175.) 

7 7 mo-v > '-O 




RECEIVER MISS DISTANCE (x - R ) (km) 


FIG. 21. REFRACTION GAIN G r (x-R ) AT 50 Me DURING 
OCCULTATION ALONG TRAJECTORY 1. 



FIG. 22. DECREASE IN PHASE PATH (s p ) DUE TO THE 
MARTIAN IONOSPHERE (f = 50 Me, TRAJECTORY l). 


SEL-64-025 


42 



FIG. 23. RATE OF CHANGE IN PHASE-PATH DECREASE 
dSp/dx FOR TRAJECTORY 1. 


the Martian ionosphere causes a doppler shift of 1.7 cps at miss distance 
600 km. As explained in Sec. IIB, this can be measured by comparing it 
with the eighth subharmonic of a 400-Mc signal. 

The shape of the curves for phase path and refraction gain vs re- 
ceiver miss distance depends strongly on the trajectory. Figures 24 and 
25 show the results for trajectory 2. We see that the gain and the phase 

path now are multivalued functions of (x - R ) . This is because the 

P 

signal in certain regions is received over more than one propagation 
path at a time. 

Trajectory 2 crosses two branches of the outer caustic at miss 
distances 1700 km and 1300 km, respectively. The inner caustic in this 
numeric example has only one branch in the x-y plane. This branch is 
encountered at miss distance -1100 km. 

So far we have neglected the nonspherical variations in the electron 
density distribution. Taking these into account we find the changes 
shown in Figs. 26 and 27. The nonspherical ionosphere used here corre- 
sponds to Eq.(4.l) for X - 85 deg and 

x( — - — > y) = x( — * 85 deg) for Y > 85 deg 

\cos a V 'cos a 1 a & 
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[A modification of Eq . (4.l) is necessary because it is not valid for 
large X*1 The two ionospheres used in the comparison in Figs. 26 and 
27 have the same electron density only along the x axis. 
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FIG. 25. REFRACTION GAIN G (x-R ) FOR TRAJECTORY 2. 

r p 



RECEIVER MISS DISTANCE (km) 


FIG. 26. EFFECT OF NONSPHERICAL PERTURBATIONS IN 
THE ELECTRON DENSITY DISTRIBUTION (s AT 50 Me 
FOR TRAJECTORY l). P 
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FIG. 27. EFFECT OF NONSPHERICAL PERTURBATIONS 
IN THE ELECTRON DENSITY DISTRIBUTION (REFRAC- 
TION GAIN AT 50 Me FOR OCCULTATION ALONG 
TRAJECTORY l). 
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V. DETERMINATION OF THE ELECTRON DENSITY DISTRIBUTION FROM AMPLITUDE, 
PHASE PATH, OR GROUP PATH 

Previous chapters have discussed how to find amplitude, phase path, 
and group path of high-frequency radio waves propagating through a 
planetary ionosphere with known electron density profile. The purpose 
of an experiment would of course be to measure the gain G^(x), the 
changes in the phase path s^(x), or the changes in the group path 
8 (x) as the transmitter is occulted by the ionosphere, and from this, 
deduce the distribution of electron density for that part of the iono- 
sphere which is probed by the waves reaching the receiver. In other 
words the problem has been reversed, but the results derived previously 
are still of great value. 

Propagation through a nonspherical ionosphere can, as has been 

shown, be treated as if the ionosphere were spherically symmetric. We 

then use the same radial distribution of electron density in that part 

of the ionosphere which is probed by the signal, but variations in the 

electron density with longitude and latitude are neglected. 

The equivalent trajectory (Fig. l) comes in very handy here too; 

and we will assume that G (x), s (x), and s (x) correspond to 

r p g 

measurements along this trajectory. 

The two first sections in this chapter are mainly devoted to 

studies of processes which can be used to find the angle of refraction 

2t(r ) or the normalized straight-line integrated electron density 
o 

s (r ) from G (x), s (x) or s (x) . Here x denotes the abscissa 
2 N o ' r ' p 7 g 

coordinate of the receiver and r the radius of closest approach for 

o 

the raypath. (See Fig. 5.) 

Later we show how the normalized electron density distribution 
X(r) can be found from s 2 (r o ). 

A. DETERMINATION OF THE NORMALIZED STRAIGHT-LINE INTEGRATED ELECTRON 

DENSITY s 0 (r ) FROM THE REFRACTION GAIN G (x) 

& o r 

The refraction gain G (x) along the equivalent trajectory is given 
by Eq. (2.33). We will assume that that part of the equivalent tra- 
jectory along which the occult at ion is observed can be approximated with 
a straight line: 

y = y^ - UX (5.1) 
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The same procedure can be employed when a piecewise linear approximation 
is used for the equivalent trajectory. 

The radius of closest approach r^ for the ray that is crossing the 
receiver will be approximated by 

r = x - (2t + P) y (5.2) 

O 

where (x,y) are the receiver coordinates. This corresponds to setting 

r equal to (R + q) in Fig. 5. We will further assume that y « D. 
o p 

The geometrical quantities used here are shown in Fig. 5. 

Some remarks are also necessary on the notation we are going to use. 

A ray with radius of closest approach, r^, is bent an angle 2^(r^) 
which is a unique function of r^. Since the rays in the x-y plane cross 
the equivalent trajectory, the same angle can also be considered a 
function of receiver abscissa x and we will then write 2^ (x). Al- 

■Xr 

though 2^(r o ) and 2f (x) are equal when they refer to the same ray, 
it is still convenient to have some difference in the notation, espe- 
cially since 2t (x) may he a multivalued function of x. For s and 

-K- -¥r ^ 

p, we will similarly use s (r ), s (x), p(r ), and (3 (x). For 

Ct O O 

6(r ) we have 
K o' 


r 



From Eq . (2.33) one can now obtain 


d[2T|f*(x) + p*(x)] 

1 + u[2t*(x) + p*(x)] 


g(x) dx 


where 


g(x) 



1 


G r (x)/10 

10 

_|_ — 

/ y - ux v / „ 

(1 + -)(l - 24T (x) 



(5.4) 


(5.5) 


The upper sign in Eq . (5.5) should be used when there is no caustic 
between the receiver and the planetary ionosphere; and the lower sign 
should be used when the rays have crossed each other somewhere between 
the receiver and the ionosphere. The expression for g(x) can be 
simplified when the trajectory lies close behind the ionosphere. 


SEL-64-025 


48 



Rewriting Eq.(5.4) in integral form and neglecting u[2\|r*(x) + p*(x)], 
one obtains 


2^ (x) + p (x) - P (x fe ) = J g(x) dx 

x b 


(5.6) 


where x^ corresponds to a position of the receiver before the occulta- 

tion has started. At x = we have 

b 

G (x. ) = 0 

r b 

2t*(xJ = 0 

b 


• DTT--' ux. 


(5.7) 


Equation (5.6) may at first appear to be useless because g(x) 
depends on the unknown function 2\|f (x). However the term containing 
24 |r (x) in the expression for g(x) is small compared to 1 in the high- 
frequency case. This makes Eq . (5.6) a very useful tool for obtaining 
2 1 (x). 

Let us first look at a method that can be used to integrate Eq . (5.6) 

directly. The integration can be done by summing -g(x) Ax from 

large toward smaller values of x. Let us assume that 2i|/ ^x) and 
* 

P (x) are found in the region from x to x (where x is less than 

^ -x- 

x ). The next step consists of determining (x - Ax) and p (x - Ax). 

® * * 

One can calculate g(x - Ax) using 2 \|r (x) instead of 2 \Jr (x - Ax). 

* 

This is possible because the term containing 2\|r (x) in Eq . (5.5) is 

-x- 

small compared to 1. Equation (5.6) now gives [ 2xjr (x - Ax) + p 

(x - Ax)], which again can be used to calculate the corresponding r 

* ^ 
by means of Eq . (5.2). Thus P (x - Ax) can now be computed from 

* 

Eq . (5.3). This also enables us to find 2^ (x - Ax) since we already 

know [ 2 ^ (x - Ax) + p (x - Ax)]. By repeating this procedure one 
* * 

obtains 2^ (x) and p (x) for all values of x where G_^( x) is 

known. This method is well suited for use on a digital computer when 

G (x) is known in sufficient detail, 
r 
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Another way to approach this problem is to neglect 2t (x) in g(x) 
as a first approximation. One then obtains 


2t 1 *(x) + p i *(x) 


- p*«s) 



g ] _(x) dx 


(5.8) 


where 


e 1 (x) 



1 T 


10 


G r (x)/10 


1 + 


y o - ux 
D 


(5.9) 


By using Eqs . (5.8), (5.2), and (5.3) one can now determine a first 
approximation [ 2 'Jr. (x)] for the angle of refraction 2 \| f (x). 

.y. 

A second approximation [ 2^ 2 (x)] for 2\|f (x) can be found by 
repeating this procedure. One obtains 


2+2 (x) + p 2 (x) 


p <S> ■ 


X 

-f 

x b 


g 2 (x) dx 


(5.10) 


where 


g 2 (x) 


V - ux 
“ o 


1 T 


10 


G r (x)/10 


/ y - ux v r » y - ux 

(i + -2-g )(l - 2^ (x) -2-j 


(5.11) 


This iterative process can be employed until sufficient accuracy is 
obtained for 2t (x) and (3 (x). 

y 

The functions G (x), g(x) and [2^ (x) + f3 (x)] are multivalued 
r 

functions of x when the rays cross each other somewhere between the 
receiver and the planetary ionosphere. Figure 28 illustrates how 

y _y. y 

[2^ (x) + (3 (x) - (3 (x^)] is obtained from integration along the 

branches of g(x). The poles in g(x) coincide with the poles in 

G (x), and thev lie where the trajectory crosses caustics. However, the 
r 

amplitude never goes to infinity at caustics such as G (x) predicts. 

r 

(See Appendix A.) Thus the true gain G^( x) does not have any poles. 
One can therefore not expect to obtain an accurate determination of 
[2v|r (x) + 3 (x) - p (x )] by substituting G (x) for G r (x) when the 
trajectory is crossing one or more caustics. However if the detailed 
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FIG. 28. INTEGRATION OF EQ. (5.6). 


behavior of G t (x) is known around the caustics, one may from this 
estimate the corresponding ray theory gain G^(x) and thereby improve 
the accuracy. (Multipath propagation can be avoided by choosing the 
frequency high enough such that the caustics will lie behind the space- 
craft trajectory. Estimates can be made by using the results derived 
in Chapters III and IV.) 

The next step in the process of obtaining s^r^) froin G r ( x ) 

consists of finding 2\|r(r Q ). This is done by changing the variable 

from x to r using Eq . (5.2). 
o 

Finally, Eq. (2.28) relates 2\|r(r o ) and s.fr^. Thus, 

00 

s (r ) = — f 2\|r(r ) dr (5.12) 

2 o' c I o' o 
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We have now seen how s ri (r ) can be deduced from the refraction 

2 o 

gain G (x). The function s (r ) is only related to that part of the 
r z o 

ionosphere through which the signal is propagating; and we have there- 
fore eliminated the dependence upon the particular trajectory along 
which the measurements happen to be made. 

The next section shows that s (x) or s (x) also can be used to 

P g 

obtain s^r^) in the high-frequency case. 

The radial distribution in electron density x(r) in that part of 
the ionosphere which is probed by the signal, can be determined from 
s (r ). (See Secs. C and D. ) 

2 O 

Losses due to electron collisions cause absorption at lower fre- 
quencies. This absorption coefficient is proportional to the product 
of collision frequency and X [Ref. 8]. One may, by observing the 
amplitude at different frequencies, be able to determine the radial 
variation in the collision frequency (if the data are good enough). 

The highest frequency would yield X(r) and this enables us to compute 
the refraction gain at the lowest frequency. The difference between 
the observed gain and the computed refraction gain at the lowest fre- 
quency yields the radial variations in the straight-line integrated 
absorption coefficient. The problem of finding the radial variation in 
collision frequency is analogous to determining X(r) from s (r ). 

2 O 

Sections C and D can therefore be used in both cases. 


B. 


DETERMINATION OF THE NORMALIZED STRAIGHT-LINE INTEGRATED ELECTRON 

DENSITY s_( r ) FROM THE PHASE PATH s (x) OR THE GROUP PATH s (x) 
2 V o' p g 


The only difference between the form of S p (x) 


difference in sign in front of s 


1 ‘ 


and s (x) is a 
g 

[See Eqs . (2.20)and (2.23).] We 


shall first look at a method in which this sign difference is of no 


importance as far as the determination of s (r ) goes. We will use 

2 O 

s (x) here but the method also applies to s (x). 

p g 

Equation (2.20) gives 


-s^x) + s 2 *(x) = s p (x) (5.13) 

where the again is used to emphasize that s here is considered a 

2 

function of the receiver abscissa x. One can neglect s^(x) when the 
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receiver is flying close to the planetary ionosphere. Thus as a first 

approximation s 0 _(r ) to s 0 (r ), one obtains: 

21 o 2 o 

S 21 (r o } = S p (r o ) (5 ' 14) 

where the receiver abscissa x is simply replaced by the radius of 

closest approach r for the ray that crosses the receiver. The first 

approximation in 2 ^ f 1 ( r c> ) can now be determined from Eq . (2.28). 

The function s^,(r ) may not be a sufficiently accurate approxi- 
2 1 o 

mat ion to s^r^), but this a PP roxiinat ion can be improved by repeated 
use of Eq . (5.13). For instance, when the waves can be considered plane 
before they are refracted in the planetary ionosphere, this gives: 


S 22 (r o ) = S p (x > * c <X ' ‘'o ) 




(5.15) 


This equation is obtained by combining Eqs . (5.13) and (2.24). The term 

in Eq . (2.24) has been neglected and so has the difference between 

a and r . When the receiver trajectory is given by Eq . (5.l), one can 
o 

show that 


x 


2 t- i (r o ) 
1 + u 2^i"(r o ) 


r + 
o 


This enables us to determine the second approximation S 22^ r o^ of 

s (r ) from Eq . (5.15). 
z o 

Higher order corrections to s (r ) can be found by repeating this 

^ o 

procedure . 

Another way to approach the problem is to consider the slope of 

q dsp 

s (x). For I 2u - ” — ~ I « 1, one can show that 

p ' I f dx 1 


2t (x) 


c 

f 


ds ( x) 
P 

dx 


(5.16) 


where 24 f (x) denotes the angle of refraction as a function of receiver 
abscissa x. This equation can be considered a generalization of 
Eq . (2.28). Once (x) is determined from Eq . (5.16), one can use 

the method explained in Sec. VA to find s 2^ r 0 ^‘ [ Ec l uation (5.16) does 
not apply to s^(x).] 
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Other ways to determine s (r ) have been found during this study. 

^ o 

However, they are believed to be of less importance and are therefore 
not given here. 

We have shown that either the phase path, the group path, or the 
amplitude can be used to obtain s^r^). However, it may be an advantage 
to make two of these measurements simultaneously in order to reduce 
experimental errors and possible uncertainties introduced by variations 
in the medium between the transmitter and the planetary ionosphere. 
Ambiguities can also be avoided by receiving on several frequencies in 
t he spacecraft . 


C. DERIVATION OF A RECURSIVE FORMULA RELATING THE NORMALIZED ELECTRON 
DENSITY DISTRIBUTION X(r) TO THE NORMALIZED STRAIGHT-LINE INTE- 
GRATED ELECTRON DENSITY s ( r ) 

2 o 

It was shown in the two previous sections how G^(x), s (x), or 

s (x) can be used to obtain s rt (r ). In this and the next section we 
g 2' o' 

will show how to determine the normalized electron density distribution 

X( r ) from s ( r ) . 

z o 

In this section the coordinate system shown in Fig. 29 is used. 

Let us consider the planetary ionosphere built up of K spherical 
layers with constant electron density in each. Equation (2.27) can now 
be rewritten in the following form: 

m-1 

2 j Sg( m) = X(m) Ay(m,m) + 2^^ X(n) Ay(m,n) (5.17) 

n=l 


Here s is denoted s (m) for the m^ ^ ray and X(n) is the 
^ ^ t h 

ized electron density in the n layer. In Fig. 29, Ay(m,n) 
nates the line element corresponding to the n^* 1 layer and the 
From Eq . (5.17) , 


X(l) 




Ay(i, 


(i) 

FT 


normal - 

desig- 

th 

m ray . 


X(2) 




1 



Ay(2, 


X(n) Ay( 2 , n) 

2 ) 
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(5.18) 


x(m) 


2 f s 2 (m) 


m-1 

- 2^ X(n) Ay(m,n) 
Ay(m,m) 


which provides us with a recursive formula relating X to s 

The two-dimensional matrix Ay can be expressed in terms of the 
radii of the spherical layers: 


Ay(m,n) = 
for m > n, and 


2 / r m + r m4l \ 2 

r n ' V 2 / 


n 1 / 2 


21 1 / 2 


2 . r + r _ 

r , + / m m+1 
n+1 


Ay(m,m) = 2 


m 


r + r 
m m+1 


1/2 


(5.19) 


(5.20) 


for n = m. 


It turns out that these results can be extended to include a special 
case where the ionosphere is not spherically symmetric. Although this 
case is not of much importance in this context, it may clarify the 
inversion process and will therefore be included. 

Let X be a separable function of the radial distance to the center 
of the planet r and of the solar zenith angle X : 


x(r,X) = X x (r) A(X) 


Using A(m,n) and A'(m,n) to denote the solar zenith angles for the 
line elements Ay(m,n) above and below the x axis respectively, one 
finds 


X 1 (m) 


m^l 

2 | s 2 (m) X 1 (n) Ay( m, n)[ A(m, n) + A'(m,n)] 

lfel 

Ay( m,m) A(m,m) 


(5.21) 


which shows that it is necessary to know A(x) in order to find X^(r). 

Another convenient method for calculation of X(r) from S 2^ r 0 ^ 
is given in the next section. 
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FIG. 29. IONOSPHERE CONSISTING OF K 
SPHERICAL LAYERS. 
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D. APPLICATION OF THE ABELIAN INTEGRAL EQUATION TO THE DETERMINATION 

OF THE NORMALIZED ELECTRON DENSITY DISTRIBUTION X(r) FROM THE 

NORMALIZED STRAIGHT-LINE INTEGRATED ELECTRON DENSITY s„(r ) 

2 o ' 

Equation (2.27) is a special case of an integral equation solved by 
Abel. His results can be found in most textbooks on the subject [Ref. 9]. 
Changing the variable of integration in Eq . (2.27) gives 


c 

f 


s (r ) 
2 V o / 


00 




(5.22) 


2 

where r is given by Eq. (2.9). By using l/r as a dummy variable 
of integration in Eq . (5.22), one can obtain from comparison with Abel's 
solution: 


X(r) 


2 c _d_ 
ji f dr 



(5.23) 


This is the formal solution to the problem, but it is not a very 
convenient solution from the point of view of numerical analysis. It 
is therefore tempting to try to modify Eq . (5.23). One can show, after 
differentiation with respect to r, that: 



2 2 " 3 / 2 

s (r )] r (r - r ) dr 

2 o' o o ' c 


(5.24) 


This form of the integrand cannot be evaluated accurately close to 
the lower limit r. Expanding the integrand in the interval r to 
r + Ar gives : 


X(r) =- 


c dS 2 (r) 


dr 


( 

\ 2l \ 


1/2 


2 c f 

n f J 

r+Ar 


U 2 (r) - 


s 2 (r o )] 


-Q A> 


r - r 


r dr 
o o 


(5.25) 
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where Ar should be small enough so that s^r^) can be a PP roxima ted 
accurately by s 2 (r) + [ds 2 ( r ) /dr ] ( r^ - r) . It is also assumed that 
Ar « r . 

The function s (r ) in the integrand of Eq . (5.25) decreases 
^ ° 2 2~3/2 

rapidly with increasing r , but r (r - r ) does not . There- 

o 0 0 

fore, in order to cut down the computation time, one can modify the 
integral in the following manner. Let 

s rt (r ) ~ 0 for r > r, 

2 o o b 


Splitting up the integration in two new intervals (r + Ar, 
(r , 00 ) gives 


X(r) 


2 c 

r. f 


ds 2 (r) 

dr 



r+Ar 


S 2 (r o )] 



and 


•( 


2 

r 

o 



r dr 
o o 


2 c 
K f 




(5.26) 


and this form is very well suited for use on a digital computer. 

We have now discussed procedures which can be used to obtain the 

normalized electron density profile X(r) from measurements of changes 

in the amplitude, the phase path, or the group path during occult at ion. 

Because of the complexity of the methods and the many approximations 

involved, it appears almost impossible to carry out any general analysis 

of the final errors in X(r), However, let us briefly look at the 

simplest case where refraction is negligible so that s 2 (r^) Can 

determined directly from phase-or group-path measurements. We will use 

s 2E ( r o ) to denote the error in s (r ), and X g(r) denote tIie 

corresponding error in x(r). One can show that X (r) is given by 

Eq . (5.26) when s 0 (r ) is replaced by s OTr ,(r ) [and x(r), by X (r)]. 

z o Zh o & 

[Errors in s_(r ) for r > r do not influence the final result when 
2 o' o b 

Eq . (5.26) is used in the inversion process.] Equation (5.26) can be 

simplified when s (r ) can be considered as varying rapidly around 
2E o 

zero so that 
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b 

2 c f 

J 


-3/2 


s (r )(r 2 - r 2 ) r dr 

2E o /N o / 00 


r+Ar 


may be neglected. One then obtains for the absolute error in x(r); 


- 1 1 


r ds 0 J r ) 


2E 


dr 


Ar + s 2E (r) 


(2rAr) 


- 1/2 


showing that the error in X at r in this case only depends on 
s (r) around r. In order to determine the standard deviation in 
X (r), it would be necessary to assume a probability distribution for 
s 2E ( r ) and ds 2E (r)/dr. 

The experimental errors may be the largest errors, but we also like 
to check the approximations we have developed. A complete numerical 
example is therefore worked out in the next section. This example may 
also clarify some of the previous derivations. 


E. NUMERIC EXAMPLE 

A numeric example was given in Sec. IVD, where it was shown how 
amplitude and phase may vary during an occult at ion. In particular, 

Fig. 21 shows G^(x - R^) for trajectory 1. 

We will now use, as an example, the gain shown in Fig. 21 and from 
this try to determine the corresponding electron density distribution 
in the Martian ionosphere. In doing so, we will make use of the results 
developed in Secs. VA and VD. 

Equation (5.9) and Fig. 21 can be used to determine g^(x) . We 

note that 1 /D and p in Eqs. (5.9) and (5.8) are both zero in this 

example because in Sec. IVD we neglected the curvature of the wavefronts 

before refraction in the Martian ionosphere. (The quantities D and 

(3 are defined in Fig. 5.) The result of computing g^(x) is shown in 

Fig. 30 where x - R is the receiver miss distance as viewed from the 

P 

earth . 

Integration of Eq . (5.8), where g (x) is taken from Fig. 30, gives 

* 

2^ (x) as shown in Fig. 31. 


59 


SEL-64-025 





SEL-64-025 


60 


This procedure can now be repeated. By using Eqs . (5.1l) and (5.10), 
we obtain g 2 (x) and (x). The result is shown in Figs. 30 and 31. 

-¥r 

Neglecting the difference between 2f (x) and 2\|f (x), we can now 
determine 2^(r Q ) by means of Eq. (5.2). The result of this nonlinear 
scale transformation is shown in Fig. 32. 



FIG. 32. THE FUNCTION 2 \| r (r -R ) AT 50 Me FOR A 

MARTIAN IONOSPHERE. ° P 


The next step consists of determination of s (r ). By using Eq . 

z o 

(5.12), one obtains the result shown in Fig. 33. 

Finally, Eq . (5.26) gives the electron density distribution that is 
causing the gain fluctuations in Fig. 21. The calculated points are 
plotted in Fig. 34, together with a curve representing the Chapman model 
originally assumed in Sec. IVD (Fig. 19). The agreement between the 
assumed and the calculated electron density distribution is good in 
spite of the many approximations involved. 

Looking at both Secs. IVD and VE, we see that we have ended where we 
started, namely with the electron density distribution. When the data 
from an occult at ion experiment become available, it may be useful to 
check the result with a similar loop: 

1. Use the data to determine the electron density distribution, then 

2. Use the electron density distribution and calculate the quantities 
that were measured during the occult at ion. 
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Note that the last step can easily be carried beyond the first-order 
approximation (high-frequency approximation), while the first step 
cannot . 


s 2< r o - V 



FIG. 33. THE FUNCTION s (r -R ) AT 50 Me FOR A 
MARTIAN IONOSPHERE. 2 p 



FIG. 34. COMPARISON BETWEEN ASSUMED AND 
CALCULATED ELECTRON DENSITY PROFILES. 
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F . OTHER ALTERNAT I VES 


In the previous sections we have seen some methods that can be used 
to determine the electron density distribution from measurement of 
amplitude, phase path, or group path. However, the methods are 
laborious, and it may turn out that the data, when available really 
will not qualify for this type of procedure. In other words, it may 
be that the uncertainties will be so large that, even after careful 
smoothing, the data may yield an unrealistic electron density distribu- 
tion. In this case it may be better to determine characteristic 
properties of the ionosphere by comparing the data with the results 
obtained in Chapter III or Chapter IV. From this it may be possible, 
for instance, to determine upper and lower limits for the parameters of 
a Chapman-model ionosphere that will yield amplitude, phase path, or 
group path in agreement with the measurements. 

Faraday rotation is another quantity that can be measured as 
mentioned in the introduction to Chapter II. One can show, when the 
quasi-longitudinal approximation applies, that the Faraday rotation is 
proportional to the integral along the raypath of X times the longi- 
tudinal component of the magnetic field [Ref. 10], This yields a means 
of exploring the planetary magnetic field. However, assuming that the 
magnetic field around the planet may be considered a dipole field, one 
may also use the Faraday rotation to check or improve x(r) found from 
phase or amplitude. The dipole field is uniquely specified by three 
numbers (two angles and the dipole moment), and the rest of the informa- 
tion contained in the Faraday rotation vs miss distance curve is there- 
fore on x(r). 

So far we have studied radio-wave propagation through planetary 
atmospheres and ionospheres. It is also necessary to consider the 
waves reflected from the surface of the planet. The next chapter will 
deal with this problem. 
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VI. REFLECTION OF RADIO WAVES FROM PLANETARY SURFACES 


Section IVD shows how amplitude and phase may vary for a signal 
propagating through a Martian ionosphere; however, waves will also be 
reflected from the surface of Mars. Let us therefore, in order to 
complete the example in Sec. IVD, briefly consider the reflected signal. 
Figure 35 shows the ray picture around the limb of Mars, assuming a 
smooth surface. It is clear from this picture that the intensity of 
the reflected signal will be way below the intensity of the direct 
propagating signal both along trajectory 1 and along trajectory 2. 

(See Fig. 18.) It therefore appears easy to distinguish the direct 
traveling wave from the wave reflected from the surface. This pre- 
liminary conclusion is based on the simple model used in Fig. 35. How- 
ever, a more detailed study in the next sections will show essentially 
the same result. 



FIG. 35. REFLECTION FROM A SMOOTH, PERFECTLY REFLECTING PLANETARY 
SURFACE. 
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A. APPLICATION OF HUYGENS -KIRCHHOFF PRINCIPLE 


The balance of this chapter considers radio waves transmitted from 
the earth and reflected from a rough planetary surface. A receiver is 
flying by the planet* The planetary surface is moving relative to the 
receiver, and this motion will impose a phase modulation on the re- 
flected signal. 

Many authors have studied radio-wave reflections from the moon and 
planets for the monostatic case in which both transmitter and receiver 
are situated on the earth [Refs. 11, 12, 13], Here we will extend some 
of this work for use in the case where the transmitter is on the earth 
and the receiver is in a spacecraft flying by the planet (bistatic 
radar). The geometry at time t is illustrated in Fig. 36. The 



WAVES FROM 
THE EARTH 


FIG. 36. GEOMETRY OF THE PROPAGATION PATH AT 
TIME t. 
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P r ] t -coordinate system is chosen with the origin in the center of the 
first Fresnel zone and both the | and the T]^ axes are tangential 
to the mean planetary surface. The ^rj^^-coordinate system is moving 
relative to the surface such that its origin always remains in the 
center of the first Fresnel zone. We also need a coordinate system 
that is fixed to the surface. Let us call its axes T] , and £. 

The last coordinate system coincides with the ^r^f^-coordinate system 
at time t as shown in Fig. 36. 

We will apply the Huygens -Kirchhoff principle in order to try to 
relate the statistics of the reflected signal to the statistics of the 
planetary surface. In order to do that, we need an expression for the 
phase of the signal reflected from a point on the planetary surface. 
Figure 37 shows a point (£ q , T)^, £ ) on the mean planetary surface. 

The reflection takes place on the surface at a distance h(£^, T) , O 
below the mean surface. It is necessary to make approximations in the 
expression for the phase path. Assuming that h(£ Q , t ) is sma11 

compared to the other dimensions involved, one obtains for the phase 
path between the transmitter, the reflection point, and the receiver, 


k(d + d)«k 
r t 


2/2 2 , 
£ / cos 9 cos 9 2 cos 9 N 

...of O o o 

d + d + _ 1 + — - + — 

or ot 2 \ d d , R , 


or 


ot 


T1 / . - 2 cos 9 \ 

+ (-T — + -z — + -) + 2h cos 9 

2 \d d ^ R / o 

\ or ot p / 


( 6 . 1 ) 



MEAN PLANETARY 
SURFACE 


FIG. 37. CROSS SECTION OF SURFACE FEATURE. 
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where k is the free-space wave number, is the mean planetary 

radius, G is the angle of incidence at the origin, and d and 
o or 

d^ are the distances between the origin, the receiver, and the trans- 
mitter, respectively. Only the first-order terms are maintained in 
Eq . (6.l). We have also neglected the effect of the ionosphere because 
this only causes a constant phase shift when the ionospheric contribu- 
tion is independent of the position on the surface where the reflection 
takes place. 

Let be the amplitude at unit distance of the electric field 

of the transmitter. The Huygens -Kirchhoff principle then gives for the 
amplitude of the reflected wave at the receiver: 


E (t) = 
r 


£//■ 


dS 

exp { i [ ojt - k(d + d )ll(cos 9^ + cos G ) - — — 
t r 7 ' t r' d d 

t r 


(6.2) 


where 9^ and 6 are the angles that the normal to the surface 
t r 

element dS makes with d^ and d^ , respectively. The Huygens - 
Kirchhoff principle applies when the surface does not have any roughness 
scale comparable to the wavelength. 

Equation (6.2) can be simplified by neglecting the variation in the 
distance and obliquity factors over the region of the surface con- 
tributing to the reflected signal. This simplification gives 

-+00 40c 


E r (0 


ikE cos 0 r r 

2nd ff 

or ot J J 


exp {i[ut - k(d + d^)]} d| Q dr} Q 


(6.3) 


where the limits in the integration have been set equal to ± 00 , This 
is valid when a large number of complete Fresnel zones take part in 
the reradiation from the surface. (This approximation is not valid 
when the first Fresnel zones approach the shadowed part of the 
surface . ) 

We will define the time autocorrelation function r(t) of E (t) 
in terms of an ensemble average: 


R( 0 = 6/(0 E r (t + T)) 

o ' 


(6.4) 
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where z is the free-space wave impedance. The dimension of r(t) 
o 

is power per unit area. 

We can now combine Eqs . (6.3) and (6.4), which gives 


r(t) 


kE^ cos 0 

2jtd cT7 
or ot 


+00 +oo 


*)< // 


exp { i [ -ut + k(d t + d r )]) d| Q dri Q 


+00 +00 

// 


exp 


{i[u)t + ujT - 


k( d 1 + d 1 ) ] } d£ 1 
\ t r o 


dT] 1 
o 


(6.5) 


Note here that the last double integral is not equal to the first double 
integral since the coordinate system used for the integration has moved 
during the time T. 

For the phase path in the last integral one can show that 


k(d ' + d *) 
t r 


2 2 2 
£ 1 /cos 0 cos 0 

/ • • v o o o 

d + d + ( d + d ) T + — — I — + 

ot or ot or 2 \ d d 


or 


ot 


2 cos 0 \ q 
o \ o 
+ + 


R 


+ 2h 1 cos 0 


„ , 2 cos 0 ' 

11 o 

+ - — + 


2 l d d 
\ or ot 


R 


( 6 . 6 ) 


where d and d denote the time rate of change of d and d , 
ot or ot or 

respect ively . 

Combining Eq . (6.5) and Eq . (6.6), one obtains: 


2 +00 +50 +00 +00 

/kE COS 0 \ § 

R(T) d~° exp ( i[wT ' k( V + d ot )T]) 

o \ or ot / 

— ocT —co —oo —oo 


ffff 
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2 
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(6.7) 


exp 


i 

1 2 ^ d 


or 


ot 


2 cos 0 c 
R 


/ 2 f 2 v 

(lo ‘ V ) 


-/exp [12k cos e o (h - h’)]) d! 0 dS 0 ' dn o di) o ’ 


where 


h = h(s o , n c , t) 


h’= h(l o ', T| o ', t + t) 


The next step in this development consists of changing dummy 
variables of integration. To be more specific, what we want to do is to 
integrate in a coordinate system that is fixed to the surface instead 
of moving along with the Fresnel zone pattern, such as our present 
coordinate system does. We have already defined a useful coordinate 
system for this purpose in Fig. 36. This was the |r]£-coordinate system 
which coincides with the subscripted coordinate system (! o *n o £ system) 
at time t (see Fig. 36), but is fixed on the surface. 


Let 


and 


be the velocity with which the center of the 


first Fresnel zone is moving over the surface in the £ and r\ 
directions, respectively. Then we have at time t, 


6 = t. 


n = ti 


and at time ( t + t) , 


5* = V T + £ q ' 


ti 1 = V T + r> 1 
T) l O 


If we also introduce A£ = £ f - £ and Ar| = T| 


= n f - 




we can use £ , r] , AS, and AH as dummy variables of integration. 
This gives 

5 - s t -v i 


= n + An - v t 


( 6 . 8 ) 
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! o " ^o’ 2 = (v | T " A ^ 2 ^ + A ^ - V | T ) 
\ - n o ' 2 = (v^T - At] ) ( 2t] + ATI - V^T) 

and h - h' = h(|,ri) - h(^ + Af , n + Ari) 


(6.9) 


( 6 . 10 ) 


We will now assume that h has a gaussian joint -probability density 
function. This gives 


< \ 2 2 2 / 2 2 
exp [ i2k cos Q Q (h - h')]> = exp {- 4k h cos e^l-pWA^ + Ar\ )]) 

( 6 . 11 ) 

(The left side can be considered the joint characteristic function of 
the gaussian random variables h and h 1 . Expressions for the joint 
characteristic function of gaussian random variables can be found in 
most textbooks on the subject.) 

In Eq. (6.1l), h 2 denotes the mean square height and pWAp * +Atj 2 ) 
is the normalized autocorrelation function for h. The mean of h is 
zero since we have chosen the mean planetary surface with radius 
as the reference level. 

We can now combine Eqs . (6.7), (6.8), (6.9), and (6.1l), which 
gives : 
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+00 
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cos 9 cos 9 2 cos 9 \ 
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exp 


. 1 1 
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2 cos 0 


R 


(v^T - AT) )r) 


dr| d(A|) d(An ) 


( 6 . 12 ) 

The integrations in | and T] yield 5 functions. After inte- 
gration also with respect to and A], one finds 

i f o o o 

r(t ) ^ R(o)exp { i [ ojT - k(d + d )t]} exp { - 4k h cos 9 [l - p( vt)] } 


where 

and 


V = si + V 2, 

i 


r(o) ^ — 


cos 0 


z \ d d ^ J / cos 0 cos 0 
o \ or ot / / o o 2 

Id + d . + R 

\ or ot p 


)(^ + 3T 


2 cos 0 


(6.13) 

(6.14) 

(6.15) 


Equation (6.13) shows that by measuring R(t) we can determine 
p(vT) if the rest of the parameters are known. 

We see that h^ = 0 gives a doppler shift of the reflected signal 
but no frequency broadening. However, we are not interested in smooth 

p n 

spheres but in rough planets. For h » A one has the situation 
illustrated in Fig. 38, where it is assumed that the derivative of 
p(vT) exists for vT = 0, and that for small (vt): 


p(vT) « 1 -(£-) 


(6.16) 


This gives 


R(t) rs r(o) exp { i[ u) - k (d Qr +d ot )] T ) exp 4k 2 h 2 cos 2 0 o ^-^r-j J 


(6.17) 
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In this case we see that it really does not make any difference what 
shape p( vt) has far away from the origin as long as it is less than 1. 
The reason is that the factor (4k h 2 cos 2 ©^) has such a large value 
that r(t) for practical purposes is zero, unless p(vr) is very 
close to 1. Under these conditions we can only hope to determine p(vr) 
from r(t) for small values of the argument, and this may not be enough 
to draw conclusions about surface-roughness components with larger 
scales . 


B. SPECTRAL BROADENING OF THE REFLECTED SIGNAL 

Equation (6.17) gives, for the power spectral density of the re- 
flected signal, 


s(f) 


Alr(o) 


4 si it hr v cos 0 


exp 


_L ( cL 
^2 \ 4v cos 




where 


w 1 (a i \ 

= t r ’ =r U x d ) 
2ti A ot or ; 


(6.18) 


= center frequency in the spectrum of the reflected signal 
c = free-space phase velocity 

A - free-space wavelength at the center frequency 

L = horizontal scale parameter for the surface irregularities 

0 = angle of incidence at the center of the first Fresnel zone 

o 

v = the velocity with which the center of the first Fresnel zone 
is moving across the planetary surface . 

The power spectrum takes on other shapes for other values of h 2 . 

For h 2 equal to zero, the spectrum becomes a 5 function; and for 
intermediate values of h 2 , one can find an expansion for the spectrum 
containing a 5 function plus an infinite series [Ref. 14]. The power 
associated with the 5 function has been assumed negligible in Eq . (6.18). 

Finally we obtain the rms bandwidth of the spectrum from Eq . (6.18): 

__ v cos 0 O 

2Af = 4 \l~2 h2 f c — — (6.19) 
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The rms bandwidth of the reflected signal can also be found directly 
from Eq. (6.13) without first assuming that h 2 » A 2 , and that p(vT) 
for small arguments is given by Eq . (6.16), Assuming again that the 
initial slope of p(vr) is zero, one obtains: 


2Af 


4v cos 0 


h 2 


d 2 p(x)\ 

X ) 


11/2 


x=0 


( 6 . 20 ) 


C. THE RADAR CROSS SECTION ( cr ) OF THE PLANET 

r 

The average reflected power per unit area at the spacecraft is R(o), 

given by Eq . (6.15). Since l/d ^ is negligible compared to l/d 

' ot or 

and 2 /R , we have 
' P 


R(O) 



cos 



e 

o 


2d cos 0 \ 



( 6 . 21 ) 


We can now relate R(o) to the radar cross section cj- of the 
planet. Here the radar cross section will be defined so that cr^ 
times incident power per unit area on the planet gives a power which, 
if scattered isotropically, would give the time average power per unit 
area R(o) reflected to the receiver. Thus 


R(0) 



( 6 . 22 ) 


where r^ is the distance between the spacecraft and the planet *s 

center. Figure 39 shows how the normalized radar cross section varies 

with r /R and 0 . The limiting cases r /R =1 and r /R = oo 
s' p o ° s' p s' p 

correspond to reflection from a plane and a distant sphere, respectively. 
Figure 40 shows the same result as a function of Ot, which here is the 
angle between the transmitter and the receiver as viewed from the 
center of the planet. 

It is of interest to check the results derived here by introducing 
h 2 = 0 and comparing it with the expressions one can derive for a 
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FIG. 39. RADAR CORSS SECTION OF THE PLANET VS 

r /R AND 0 . 
s' p o 
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4.0 



FIG. 40. RADAR CROSS SECTION OF THE PLANET VS r /r AND a. 

s' p 


2 

smooth conducting sphere. We notice the a does not depend on h 

at all, at least not to the degree of approximation used here. This 

lack of dependence must mean that the cr derived here is the same as 

r 

for a smooth conducting sphere. Comparison with expressions obtained 

earlier for this case shows that this is indeed true [Ref. 15], 

We assumed that the entire surface was illuminated when we applied 

the Huygens-Kirchhof f principle. This assumption is no longer true for 

large 0^ because part of the surface will then lie in the shade of 

hills and mountains. This shadowing will tend to reduce cr^ for large 

values of 0 . 

o 
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So far we have assumed that the surface is a perfect reflector. It 
is necessary to multiply the radar cross section found above with the 
power reflection coefficient when only part of the power is reflected. 

D. NUMERIC EXAMPLE 

A numeric example may be helpful in order to understand the problem 
better. Figure 41 shows the hyperbolic trajectory for a spacecraft 
flying by Venus. 



i 

'p 


INCOMING WAVES 

FIG. 41. ASSUMED FLY-BY TRAJECTORY. 
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We shall first compare the power received from the reflected wave 
and the direct traveling wave. The power received directly from the 
earth (P ) is given by 


P = — g 
d 4* e T 


rd 


D 


(6.23) 


where P^_ = transmitted power 

= transmitter gain toward Venus 

A , = receiving aperture directed toward the earth 
rd 

D = distance between Venus and the earth. 

For the power received from Venus (P^), we have : 


P. O' A 
p t _r rr 

r - 4n g T 2 2 

D 4jtr 


(6.24) 


where A = receiving aperture toward Venus 
rr 

r = distance between the spacecraft and the center of Venus 
s 

<j = the radar cross section of Venus as viewed from the 
r 

spacecraft. A power reflection coefficient of 0.1 will be 

used for the surface on Venus. 

Assuming that A = A we find 

rr rd 


P 

r 



a 

r 


„ 2 
4jir 


s 


At positions I and II in Fig. 41, this gives 


(6.25) 


10 log 


P 

r 



-25.4 db 

(position 

-17.0 db 

(position 


We see that the total power reflected from the planet is way below the 
power in the direct signal even though the receiver is passing by close 
to the surface on Venus. 
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The spectra are computed from Eq . (6.18) and shown in Fig. 42. It 
has been assumed here that the waves from the earth have a frequency of 
50 Me in the Venus frame and "Vh^/L has been set equal to l/lO. 

It is necessary to specify some more of the parameters involved 
before we can estimate the signal-to-noise ratio for the reflected wave. 
Let 

= 400 kw 
g T = 360 
g r = 1.4 

10 

D =6*10 meters. 


These values give, for the average reflected power received at the 
spacecraft , 


_ 0.366 • 10 16 w (position i) 

r 2.594*10 ^ w (position II) 

From Eq. (6.18) one can now find the maximum power spectral density: 

— 18 

0.474*10 w/cps (position I) 

“18 

0.55*10 w/cps (position II) 


S(I C ) ~ 


The total power received is increased as the spacecraft moves from 
position I to position II, but the rms bandwidth of the spectrum is 
increased too; as a result, the maximum power spectral density 
s( f ^) is only slightly increased in this example. 

Assuming a constant noise temperature of 6000 °K, one finds for 
the noise power spectral density N 


-20 , 

N = kT = 8.275*10 w/cps 

This gives the signal-to- noise ratio: 


s(f c ) 

10 log — ~ 


7.58 db 
8.22 db 


(position i) 
(position II) 
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which is proportional to L/vh^. The above signal-to-noise ratios are 
computed assuming L/ >/h^ = 10. With L / \fh^ = 4 one would get a 
reduction in the signal-to-noise ratio by 4 db, while 20 would give 3 db 
increase in the signal level. 

The power reflection coefficient in this numeric example was set 
equal to 0.1 for both receiver position I and receiver position II. 

That the reflectivity is considered equal in the two cases is just a 
simplifying assumption, since both the angle of incidence (0 q ) and 
the surface material of the reflecting areas may be different. We will 
come back to this problem in the next section. 

As we have seen, the signal reflected from Venus contains informa- 
tion on the surface correlation function for the areas from which the 
signal is reflected. By observing the reflected signal over a large 
part of the orbit, one may determine how inhomogeneous the surface is 
and may find planes and mountain terrain. The bistatic-radar experiment 
may therefore offer a unique way of exploring the cloud-covered surface 
on Venus . 


E. THE POLARIZATION OF RADIO WAVES REFLECTED FROM PLANETARY SURFACES 
AND ITS RELATION TO THE ELECTROMAGNETIC PROPERTIES OF THE SURFACE 
MATERIAL 


We saw in the previous sections how to relate the fading of the 

reflected signal to the surface roughness. In this section we shall 

briefly look into how the polarization of the reflected wave can be 

related to the electromagnetic properties of the surface material. 

We will assume that a circularly polarized wave is transmitted from 

the earth. The plane containing the transmitter, the planetary center, 

and the receiver will be called the plane of incidence. This is the 

C £ -plane in Fig. 36. We shall denote the electric field components 
o o 

of the reflected wave, E and E . Here E . is the electric field 

ri rp n 

in the plane of incidence and E is the field perpendicular to the 
plane of incidence. 

The surface reflection coefficients are different for the two 
components of the wave. By comparing the rms value of E with the 
field of the direct traveling wave, one can estimate the reflection 
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coefficient for this wave component. We have to make use of the approxi- 
mation found for the radar cross section of a rough planet in order to 
find the absolute values of the reflection coefficients. This source of 
inaccuracy can be avoided if one instead considers the complex ratio (r) : 

R = exp (±if) (6.26) 

rp 

In general, both E^ and E are complex numbers. The sign in the 
exponent is chosen so that R is equal to the ratio between the 
reflection coefficients for the two wave components. 

We will assume that the effects of small-scale roughness and multiple 
reflections are negligible. 

The conductivity cr and the relative dielectric constant £ may 
vary over that part of the surface which is contributing to the echo. 
Assuming that the distribution in cr and £ is sufficiently narrow, 
one can show that their mean cj and are given by: 

— (T . 2^ 

£ - i *=» sin 0 

r uj£ o 

o 

where i =^Ti and £° = dielectric constant of free space. The 
relative magnetic permeability for the surface material has been set 
equal to 1. 

The size of the polarization ellipse of the reflected wave (measured 

in volts/meter) will change rapidly due to the fading. However, the 

axis ratio and the orientation of the ellipse will only change due to 

changes in 0 , <j_ and £ . 

o r 

The method indicated here permits us to map a and € . However, 
the resolution will be relatively poor unless the spacecraft is passing 
close to the surface. For multiple-layer surfaces it would be necessary 
to use several frequencies in order to determine the different parameters. 


2 

+ tan 0 
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VII. CONCLUSIONS 


In this study we have shown that |bistatic radar may be used to 
explore planetary ionospheres and surface s , j 

Transmitting radio waves from the earth and receiving them in a 
spacecraft that is moving behind a planet makes it possible to observe 
the variations in amplitude, phase path, or group path caused by the 
planetary ionosphere. It is shown that measurement of either the 
amplitude, the phase path, or the group path may be used to determine 
the radial electron density distribution in that part of the planetary 
ionosphere which is probed by the signal. However, it may be an 
advantage to do a couple of these measurements simultaneously in order 
to reduce possible uncertainties due to changes in the interplanetary 
medium or the earth f s ionosphere during the experiment. The two radial 
electron density distributions obtained from the measurements during 
immersion and emersion are again related to the atmosphere on the planet. 
This relationship will also help us to make a better atmospheric model. 
The lower atmosphere may also be studied by using so high a frequency 
that the signal is only influenced by the atmosphere. Besides being of 
scientific interest, this information will be of vital importance if it 
is decided to land vehicles on the planet. 

We have also shown that the statistical properties of the signal 
reflected from a planetary surface can be related to the surface statis- 
tics. It was shown, for instance, that the surface autocorrelation 
function can be related to the time autocorrelation function for the 
reflected signal when the surface-roughness scale is larger than the 
wavelength. This represents a generalization of results obtained 
previously for monostatic radar. We have also shown that the mean value 
for the dielectric constant and the conductivity can be found from the 
shape and orientation of the polarization ellipse for the reflected 
signal (assuming that a circularly polarized wave is transmitted). This 
method may find useful application in exploring Venus, where the surface 
is covered to a large extent by clouds. 
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APPENDIX A. WAVE-THEORY MODIFICATIONS CLOSE TO CAUSTICS AND SHADOW 
BOUNDARIES 


The ray theory is not adequate along caustics and shadow boundaries. 
In these regions the propagation can be studied by means of Huygens 1 
principle. 

The caustic is the evolute of the wavefronts that have emerged from 
the planetary ionosphere. By making use of this relationship, one can 
find a convenient power series expansion for the phase of the signal. 

When maintaining terms up to the third power, one can express the fields 
around the caustic in terms of tabulated functions. This analysis was 
first worked out in optics by Airy, but the solution can also be found 
in books on radio-wave propagation [Ref. 16]. Although the geometry may 
appear to be somewhat different, the solution takes on the same form and 
the derivation is therefore not repeated here. The analysis shows a 
finite increase in the amplitude as the caustic is crossed. 

The fields along a shadow boundary can also be determined from the 
wave theory. Here the solution can be expressed in terms of the Fresnel 
integrals, assuming that it is sufficient to maintain terms up to the 
second power in the expansion for the phase. The result is the same as 
for diffraction around an opaque straight edge (assuming that the 
reflected signal is negligible and that the straight-edge approximation 
applies), but one must take into account the fact that the ionosphere 
and the atmosphere will change the apparent position of the transmitter 
as viewed by the receiver. 

One can show that the amplitude E(x,y) as a function of receiver 
position x,y can be approximated by: 


E(x,y) 

1 /2 y ^ 00 2 

[f(5 + ?)l / •H- i [£(5 + j)( r o- I D J r?) ■ 2ns 2 (r o ) ]l 


dr 


r =R 
o p 


( A. 1 ) 
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where a factor exp [ iu)t - ik(D + y)] has been omitted. Here, s (r ) 

^ o 

is given by Eq. (2.22), k is the free-space wave number, r^ is the 

dummy variable of integration, D is the distance between the trans- 
mitter and the planet, and R^ is the radius of the planet (see Fig. 5) 

Equation (A.l) gives (l - i) for the amplitude before the occultation 

starts . 

A numeric example is shown in Fig. 43. This figure shows ampli- 
tude vs miss distance at 2300 Me (the telemetry frequency on the space- 
probe mission to Mars), 20,000 km behind Mars. Equation (2.17) is used 
for the refractive index, assuming negligible water vapor, T = 200 °K, 
and P = 20 and 0 millibars, respectively. At 2300 Me the effect 
of the ionosphere is negligible. 

Equation (A.l) can be changed into a more convenient form when the 
perturbations on the wavefronts, caused by the atmosphere or the iono- 
sphere, are small. In that case it is more economical to compute the 
difference in amplitude with and without the perturbing medium present. 
This calculation will converge faster. The amplitude without any 
perturbing medium present can be found from the Cornu's spiral. 



MISS DISTANCE (km) 


FIG. 43. AMPLITUDE VS MISS DISTANCE AT 20,000 km BEHIND MARS 
(f = 2300 Me, H = 20 km, T = 200° k). 
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APPENDIX B, DETERMINATION OF THE ELECTRON DENSITY DISTRIBUTION BASED 
ON THE WAVE THEORY 

Chapter V shows how one can determine the changes in the refractive 

index from measurement, for instance, of phase or amplitude during 

occultation. This analysis was all based on the ray theory. Equation 

(A.l) can also be inverted so as to give s (r ) (and therefore also 

^ o 

the radial variation in the refractive index). This inversion gives: 



( B . 1 ) 


for r > R , when y is independent of x. In this form the solution 
o p 

requires both amplitude and phase information. 
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